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ABSTRACT 

TeV J2032+4130 is the first extended very high energy gamma-ray source, which has remained enigmatic 
since its discovery, due to the lack of identification. We report here deep radio observations covering the 
TeV J2032+4130 field and revealing for the first time an extended and diffuse radio emission, as well as a 
remarkable population of compact radio sources. Some of these radio sources are in positional coincidence 
with X-ray and optical/IR sources. Future follow up studies of these new radio sources will likely contribute to 
solve the mystery of this extended unidentified TeV source. 

Subject headings: radio continuum: stars — radio continuum: ISM — X-rays: stars — gamma rays: observa- 
tions — 



1. INTRODUCTION 

The new generation of Cherenkov telescopes has revealed 
the existence of a population of v ery high energy gamma- 
ray s ources in the Milky Way ( Aharo nian et alj 1.2005a, 
2006). Some of these sources are well identified and are 
compact sou rces like LS 5039 (lAharonian et al.l [2005bh or 
LSI+61°303 dAlbert et alj|2006l) . whereas others still remain 
without known counterpart. Among the latter there are 
some extended TeV sources (lAharonian etalJl2002l [2005c) 
for which all attempts to find their counterparts at other wave- 
lengths have failed, and are being considered a new popu- 
lation of galactic sources. Because the lack of detection of 
a low energy counterpart it has been proposed that hadronic 
processes ins tead of leptonic o nes could be behind the TeV 
emission (e.g lTorres et alj d2004l) ). 

The so far unidentified source TeV J2032+4130, discovered 
with the stereoscopic High Energy Gamma Ray Astronomy 
(HEGRA) array of imaging Cherenkov teles copes in the di- 
rectio n of the Cygnus OB2 star association (Aha ronian et all 
2002), is perhaps the most representative within the group 
of TeV sources with un known counterpart. Oth er examples 
are HESS J1303-631 (lAharonian et al.ll2005d) . the second 
unidentified TeV source, and the unidentified sources fo und in 
the HESS galactic plane survey (Aharo nian et alj2005al) . The 
identification of any of these sources, by locating its coun- 
terpart at radio, optical or infrared wavelengths, would im- 
prove our knowledge of the emission processes at very high 
energies and would help to solve the problem of the uniden- 
tified TeV sources. TeV J2032+4130 is extended, with a ra- 
dius of 6'.2, being the center of gravity (CoG) of TeV photons 
located at a J20 oo.o = 20 /l 31 m 57?0 and feooo.o = +41°29'56"8 
with arc-minute accuracy. Its emission exhibits a hard spec- 
trum with a steady flux of about a few 10~ 12 erg cm" 2 s" 1 
(lAharonian et alJ2005dl) . The large angular size of the source, 
its steady flux and the location in the galactic plane point 
strongly to a galactic nature. The absence of any evident 
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counterpart led lAharonian et al.l d2002l) to propose as possi- 
ble origins of the TeV emission the star association Cygnus 
OB2 or a jet-driven termination shock. Hadronic interactions 
within the innermost region of the winds of O and B stars 
can also produce significant gamma-ray luminosities at TeV 
energies, compatible with those found in TeV J2032+4 130 
dTorres et alJl2004t iDomingo-Santamarfa & Torres 1 12006) . 

Despite the intensive observational campaigns undertaken 
for this source at different wavelengths, no obvious counter- 
part at radio, optical nor X-ray energies has been reported 
up to now, leavin g TeV J2032+4130 pre s ently unidentified 
dButt et al.l J20031 120061) : iMukheriee et alj d2003l) ). Here we 
report the discovery of a population of radio sources well 
within the ler radius of the TeV J2032+4130 extended emis- 
sion, presenting some of these radio sources optical/IR and 
X-ray counterparts. One or several of these objects, includ- 
ing both a extended diffuse component and compact radio 
sources, could be related to the very high energy emission. 
Their observational properties and nature are presented and 
discussed in the following sections. 

2. OBSERVATIONS AND ANALYSIS 
2.1. Radio 

Most of the data presented in this work were obtained 
through observations performed with the Giant Metrewave 
Radio Telescope (GMRT) of the National Centre for Radio 
Astrophysics (NCRA) in Khodad (India), during 2005 July 
9, 2005 September 1 and 2006 July 18. We observed at 
the 610 MHz frequency (49 cm wavelength) in spectral line 
mode with 128 channels covering a 16 MHz bandwidth. The 
calibration of amplitude and bandpass was achieved by ob- 
serving 3C 286 and 3C 48, while phase calibration was per- 
fomed through repeated scans on the nearby phase calibrator 
J2052+365. Correction of the GMRT flux densities for the 
increase in the sky temperature in the TeV J2032+4130 direc- 
tion was also taken into account. In addition, we also used 
archive data obtained during 2003 April 29 using the Very 
Large Array (VLA) of the National Radio Astronomy Obser- 
vatory (NRAO) in New Mexico (USA). The VLA was in its 
most compact D configuration and the observations were car- 
ried out at the 1 .465 and 4.885 GHz (20 and 6 cm wavelength, 
respectively). The VLA amplitude calibrator was 3C 48, and 
J2052+365 and also J2007+404 were used as phase calibra- 
tors. Both the GMRT and VLA data were processed us- 
ing standard procedures within the AIPS software package of 
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NRAO. 

The results of the radio observations are presented in Figs. 
Q]and[2] The left panel of Fig. [T]contains a wide field map of 
the TeV J2032+4130 region obtained by mosaicing four indi- 
vidual pointings with the VLA at 20 cm. Extended and dif- 
fuse radio emission is clearly seen inside the TeV J2032+4130 
lcr radius and closely surrounding the CoG of TeV photons, 
indicated by the central cross. The overall morphology of the 
main extended emission features in the field is also detected in 
GMRT maps (not shown here) retaining short spacings, con- 
firming their shape and extension. A spectral index image 
(not shown here) created using matching beam VLA maps at 
20 and 6 cm strongly suggests that there is a diffuse emis- 
sion half surrounding the CoG of the TeV J2032+4130 source 
(red dashed arc in Fig. [T]left) with non-thermal spectral index 
a ~ — 1, defined as S v otv a . We note that the diffuse emission 
seen well outside the TeV source circle is thermal, hinting to 
a differentiated origin for the non-thermal diffuse radiation 
within the TeV source circle and to a physical link between 
the diffuse radio mission and the TeV source. Another re- 
markable object in the field is a prominent double-lobe radio 
galaxy whose core is at aj2ooo.o = 20 /, 32 m 01?7 and <5j2ooo.o = 
+41°37'22", with an integrated flux density of 91 ± 1 mJy at 
the 20 cm wavelength. 

The right panel of Fig. Q] illustrates a zoomed VLA map 
at 6 cm of a region close to the TeV CoG. At this shorter 
wavelength we can clearly detect here two radio sources with 
different morphology. The northern one (VLA-N) is practi- 
cally compact, while the sourthern VLA source (VLA-S) is 
not consistent with being point-like. 

The most compact source sensitive radio view of the TeV - 
J2032+4130 field comes nevertheless from the final 610 MHz 
map combining all different GMRT observing epochs. This 
map is presented in grey scale plot in the central panel of Fig. 
|2]and provides the deepest radio image with arc-second detail 
ever obtained so far of the TeV J2032+4130 position. Ex- 
tended emission has been removed by not including visibili- 
ties corresponding to baselines shorter than 1 kA = 490 m. The 
field shown covers the entire la radius and at least six com- 
pact radio sources are detected with peak flux density above 
five times the rms noise of about 90 /iJy beam -1 . Their ob- 
served properties are listed in Table Q] 

Finally, we will also mention that a short VLA exploratory 
time proposal in A configuration was also conducted by us 
at the 3.5 cm wavelength on 2006 April 1 and 3rd. As a re- 
sult, a compact (angular size <0"3) radio source was detected 
with a flux density of 0.24 ±0.03 mJy located at aj2oooo = 
20 /! 31 m 52?686 and Snooo.o = +41°30'54"56. This corresponds 
to the small cross in the right panel of Fig. Q~](within the radio 
contours of the mentioned VLA-S source), with a Core label 
for reasons discussed below. 

2.2. Infrared 

The field of TeV J2032+4130 was observed in the near in- 
frared Ks-bwd (2.2 ^m) using the Calar Alto 3.5 m telescope 
and the OMEGA2000 camera on 2005 April 29th. Another 
Ks-bwd image of the field was taken on 2006 July 7th using 
the 4.2 m William Herschel Telescope (WHT) and the LIRIS 
imager/spectrograph. These ob servations were pr ocessed in 
a standard way as described in iMarti et ail (120061) with the 
IRAF package of the National Optical Astronomy Observato- 
ries (NOAO), including sky background subtraction and flat 
fielding. The left and right panels of Fig. [2] display part of 
these observations showing the identification of near infrared 



counterparts to some of the radio sources reported in this pa- 
per. The identification is on the basis of astrometric coinci- 
dence within a few tenths of an arc-second. 

3. RESULTS AND DISCUSSION 

An important consideration concerning the new radio 
sources reported here is their coincidence with X-r ay emitters 
in the field previously reported by other authors dButt et al.l 
I2003L 120061: [Mukherjee et al.ll2003l) . Hereafter, X-ray sources 
will be referred to using the Chandra \Butt et all ((2006) identi- 
fication numbers. We will discuss this issue for the three most 
relevant objects reported in this work, namely, the VLA-N 
source (or GMRT source #3), the VLA-S source (undetected 
by the GMRT) and the GMRT source #5. VLA-N and VLA-S 
are well detected separately only at 6 cm, but an estimate of 
their individual 20 cm flux densities is feasible by fitting two 
Gaussian components with fixed positions. 

VLA-N is well within the 1" uncertainty of Chandra #234 
X-ray source and it has also a faint Ks = 16.8 ±0.1 coun- 
terpart in our WHT image (see the right panels of Figs. Q] 
and|2]l. Combining VLA+GMRT data, the radio spectrum of 
this northern VLA source can be described by S v = (1.35 ± 
0.06) mJy[^/GHz]"°- 47±0 03 and is, therefore, a clearly non- 
thermal emitter. The fit to the Chandra #234 spectrum 
gives a hard photon index of 1.1 ±0.4, unabsorbed flux of 
~ 2 x 10~ 13 erg cm" 2 s" 1 in the range 2-10 keV, and hy- 
drogen column density ~ 2 x 10 22 cm" 2 . We note that the 
data are poor and do not allow for detailed fitting. Assuming 
the accepted distance to the Cygnus OB2 association of 1.7 
kpc, the X-ray luminosity in the range 2-10 keV would be 
~ 10 32 erg s . Similarly, the 0.1-100 GHz radio luminosity 
is estimated as L rac i = 1.0 x 10 29 erg s" 1 . These values have 
to be considered as lower limits if the source is beyond the 
Cygnus OB2 distance. 

VLA-S is also of non-thermal nature, with its radio spec- 
trum given by S v = (2.6 ±0.6) mJy[i//GHz]"°' 5±a2 from VLA 
data only. Undetected by the GMRT above 0.27 mJy (3a), 
its spectral index must be significantly inverted, i.e., with 
a > +1.0 between 49 and 20 cm. This fact is suggestive 
of synchrotron emission from an optically thick core at long 
wavelengths. Moreover, this object also appears superposed 
onto the Chandra #181 error circle (see again the right pan- 
els of Figs. Q]and[2]l. Since the number of X-ray counts is too 
low to perform a proper spectral analysis, we assume a photon 
index of 1.5, and hydrogen column density of ~ 10 22 cm" 2 , 
and the same distance as above, deriving an X-ray luminosity 
in the range 2-10 keV of about 10 31 erg s" 1 . For the opti- 
cally thin radio luminosity, we obtain L rac { = 1.6 x 10 29 erg 
s" 1 . Inside the Chandra #181 circle, there is a A"s = 1 3 .9 ± 0. 1 
near infrared source which is likely of stellar nature based 
on WHT spectroscopy showing no lines with cosmological 
redshift. We propose it to be the counterpart of the X-ray 
source, provided that the marginal Chandra #181 detection is 
real. However, the complex morphology of the VLA-S radio 
source (see right panel of Fig. [U suggests that it consists of, at 
least, two different components barely resolved in the D con- 
figuration of the array. One of them could be the Chandra X- 
ray/IR source and the other the compact radio source detected 
in our VLA exploratory observation in A configuration at 3.5 
cm. The other one is possibly related to the compact core just 
commented above, but further confirmation is yet necessary. 
The galactic or extragalactic nature of the mentioned radio 
core is impossible to establish from infrared spectroscopy due 
its extreme faintness in our WHT image (Ks = 19.2 ± 0.3). 
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Another remarkable coincidence is that of the GMRT 
source #5, whose position agrees very well with a V = 11.95 
and Ks = 9.07 ± 0.02 optical/infrared star belonging to the 
Cygnus OB2 association. This objec t is catalogued as star 
MT9 1 213 with a BOVp spectral type dMassey & Thompson I 
11991b iButt et all l2Q03h . There is also the Chandra #129 
source at the position of such a bright object (see error cir- 
cle in the left panel of Fig. With the same assump- 
tions as above, this star has an X-ray luminosity of about 
3 x 10 31 erg s" 1 in the same energy range. From its non 
VLA detection at 6 cm and excluding variability, we obtain 
a spectral index constrain of a < -0.76 that points again to a 
non-thermal emission mechanism with a radio luminosity of 
L rad (0.l- 100 GHz) = 1.7 x 10 28 erg s" 1 . Colliding winds in 
an early type binary or gyrosynchrotron radiation from mildly 
relativistic electrons trapped in the strong stellar magnetic 
field could be responsible for the observed non-thermal emis- 
sion in this peculiar spectrum star. 

We recall that all the luminosities provided here are just 
lower limits if the sources are further than Cygnus OB2. In 
fact, only the GMRT source #5 is a known member of the as- 
sociation but the others could be more distant. For instance, at 
10 kpc, these objects would present luminosities almost two 
orders of magnitude larger, likely requiring the presence of a 
compact object, and thus pointing to an X-ray binary nature 
(for a model of microquasars -X-ray bina ries with jets- power- 
ing ex tended TeV hadronic emission, see iBosch-R amon et al. 
d2005k 

At this point, we should not forget the existence of the ex- 
tended source present in the 20 cm VLA map around the TeV 
CoG. This feature covers a solid angle of about 27 arc-minute 2 
and its radius would correspond to ~ 2 pc at the Cygnus OB2 
distance. An integrated flux density of 116 ±3 mJy is ob- 
tained at 20 cm. For synchrotron emission with a = -1 at the 
Cygnus OB2 distance, its total radio luminosity between 0.1- 
100 GHz amounts to L ra d ~ 3.9 x 10 30 erg s" 1 . From equipar- 
tition consideration, one can also derive a total energy content 
of ~ 6.4 x 10 45 erg and a magnetic field of 1 .6 x 10" 5 G. The 
energetics would point to an efficient injector of non thermal 
particles, coherently with the fact that the region emits at TeV. 
The magnetic field value is consistent with a slight enhance- 
ment of the interstellar medium density (ISM) value. All this 
could be explained by a shock driven by a weakly radiative 
jet in the ISM produced by some of the point-like sources, 
perhaps possible microquasars, within the TeV region. In 
any case, the radio flux densities reported here provide an 
useful upper limit t o the low energy counte rpart of the TeV 
source that (see also Mukheriee et al. (2006) for a diffuse X- 



ray component inside the TeV source region presenting fluxes 
of ~ 10~ 13 erg cm" 2 s" 1 ) allow, and call for, further constrain- 
ing the models put forward so far. 

To summarize, our results reveal that at least three non- 
thermal radio sources, one of them being a significantly hard 
X-ray source, and a non-thermal diffuse radio component are 
detected in the line of sight towards TeV J2032+4130. These 
facts completely change the situation about this unidentified 
extended TeV source. Where no counterpart was seen before, 
now several radio sources have emerged and future observa- 
tions should decide whether or not one of them is related to the 
TeV emission. All radio/X-ray sources have significantly neg- 
ative spectral indices suggestive of optically thin synchrotron 
emission. Two of these sources are also of likely stellar origin 
and one is for sure a bright early type star in Cygnus OB2. Ex- 
cluding this last case, the distance to the rest is not yet clearly 
established. Thus, their X-ray luminosities may be compara- 
ble to those of X-ray binaries if placed beyond Cygnus OB2. 
Is one or more of these objects similar to the gamma-ray bi- 
naries LS 5039 dParedes et alJl2000l: lAharonian et alj|2005bl) 
and LS 1+61303 ([Albert et al.N2006l) ? With a farther location 
in the Galaxy, our sources would have similar radio and X- 
ray luminosities and, if powering the TeV extended radiation 
via, e.g., shock with the ISM, the energetic reservoir would 
be similar as well (see for instance Pare des et al.l d2006l) ). An 
in depth exploration of all these scenarios is beyond the scope 
of this letter and will be addressed in a future work. 
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RIGHT ASCENSION (J2000.0) RIGHT ASCENSION (J2000.0) 

FIG. 1 . — Left. Large scale VLA mosaic, at the 20 cm wavelength (1.4 GHz) in D configuration and uniform weight, of the unidentified 7-ray source TeV - 
J2032+4130. The CoG of TeV emission and its statistical error are indicated by the central cross, while the big yellow circle ilustrates the lcr radius of the TeV 
extended emission. Inside it, diffuse radio emission is clearly seen. The red half arc traces the extended emission around the TeV CoG position. Color scale goes 
logarithmically from above to 126 mJy. The synthesized beam is shown at the top right corner as an 40"7x40"2 ellipse, with position angle of 15?3. Right. 
Zoomed map of an interesting area inside the TeV extended emission circle where a complex of VLA radio sources has been detected at the 6 cm wave length (4.8 
GHz). Here, the red circles indicate the location of Chandra X-ray sources in the field, labeled according to the X-ray identification number given bv lButt et al.l 
12006). The Chandra source #234, at «j2000.0 = 20' ! 31'"5L84 and <5j2000.0 = +41°31'18"84, is a hard X-ray source and is one of the strongest sources in the 
field of TeV 2032+4130. This source has also been detected at 610 MHz (source #3 in Table 1). The Chandra source #181, aj2000.0 = 20' 1 31"'52?08, <5j2000.0 = 
+41°30'51"12, is a weak source. The small cross shows the location, OJ2000.0 = 20*31"'52f686 and <5j2ooo.o = +41°30'54"56, of a compact radio core detected 
at 3.5 cm during a follow up run with the VLA in A configuration (not shown here). Uniform weight was used to create this 6 cm map with a synthesized beam 
of 1 l"9x 1 1"8 with position angle of — 13?0. The color scale goes linearly from nearly zero to 0.37 mJy beam -1 . 




RIGHT ASCENSION (J2000.0) RIGHT ASCENSION (J2000.0) RIGHT ASCENSION (J2000.0) 

FIG. 2. — Center. The central panel shows our deep GMRT map of the TeV J2032+4130 field at the 49 cm wavelength (610 MHz). The central cross and the 
big circle have the same meaning as in Fig. [T] Up to 6 compact sources are detected above a 5a peak flux density level inside the TeV extended emission circle. 
Those with both X-ray and near infrared counterparts are shown in detail in the panels aside. Uniform weighting, excluding baselines shorter of lkA, has been 
used to produce this GMRT map with a synthesized beam of 5"0x4"8 and position angle of 33?0. The grey scale is linear and goes from to 2. 15 mJy beam -1 . 
Left. Zoomed view of GMRT source#5 coincident with a bright optical/near infrared early type star according to CAHA A".v-band observations. This object also 
has a Chandra X-ray counterpart indicated by the small red circle. Right. The same area as in the right panel of Fig. [T]containing the bright GMRT source #3. 
This object has an excellent match with one of the 6 cm VLA and Chandra sources in the subimage. The radio contours are superimposed onto a WHT A'.s-band 
image showing that there is also a faint near infrared counterpart. 
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TABLE 1 

GMRT 610 MHZ RADIO SOURCES ABOVE 5<T PEAK LEVEL. 



Id. 


"J2000.0 


<5j2000.0 


Peak flux 


Integrated flux 








density (/iJy) 


(density /iJy) 


1 


20''31"'42?90 


+41°29'42?'7 


478 ±83 


478 ±144 


2 


20*31 m 5Or07 


+41°32'08"4 


454 ±85 


501 ±155 


3 a 


20*31 m 51?78 


+41°31'18"3 


1148 ±82 


1749 ±189 


4 


20*32"' lOf 98 


+41°32'39"6 


521 ±85 


504 ±142 


5 b 


20*32 m 13f09 


+41°27'24;'5 


525 ± 82 


718±178 


6 


20*32"' 16f29 


+41°30'55"6 


566 ±82 


772 ±176 



a Coincident with VLA-N, Chandra #234 and near IR counterpart. b Coincident with Chandra #129 and stellar optical/near IR counterpart, 



